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ABSTRACT 



The vortex valve pneumatic oscillator is described. Data 
obtained for variations of inlet area, bias flow and temperature 
are presented. The input impedance of the valve is determined 
experimentally and found to have a region of negative resistance. 
If the valve is supplied through a capacitance (a volume) the 
combination produces a limit cycle oscillation where pressure in 
the volume is a function of time. 

Theory predicts the frequency of oscillation to be insensi- 
tive to uhe supply pressure. Analysis of the supply pressure 
versus frequency data shows the oscillator is nearly pressure 
insensitive . 

Theory predicts the frequency of oscillation to be directly 
proportional to temperature* The static valve characteristics 
were found to vary with temperature . Therefore , the temperature 
sensitivity of the oscillator was predicted from the static 
valve characteristics. The ratio of predicted temperature sensi- 
tivity to actual temperature sensitivity was 3*2 to 1. 
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CHAPTER I 

INTRODUCTION AND OBJECTIVES 



1.1 Introduction 

The vortex valve is a pure fluid modulator or amplifier. Tho 
operation of the device depends on the dynamic interactions between 
streams of fluid rather than moving parts . The absence of moving 
parts and the absence of interface problems make the pneumatic 
oscillator attractive . 

However, the sacrifice made for having no moving mechanical 
parts is difficulty in designing the geometric shapes analytic- 
ally. In fact, an analytical procedure for complete design does 
not exist . Therefore, a great deal of the design is accomplished 
emperically. 

1.2 Objective 

The objective of this thesis is to design, construct and test 
a pneumatic oscillator whose frequency of oscillation is insensi- 
tive to supply pressure but sensitive to air temperature . Design- 
ing a device which is pressure insensitive eliminates tho need for 
regulation of the supply pressure . A fluidic temperature sensor 
is highly desirable for temperature sensing in systems operating 
with oxidizing gases or gases at extremely 'nigh temperatures. 



a component 
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An ideal application for such a device would be as 
of a fluidic speed control system for a gas turbine, 

Tho plan of attack was straightforward: design an oscillator 

that is theoretically pressure insensitive and temperature sensitive , 
vary one geometric parameter to provide more information for future 
designs and test the device. 



1.3 The I-ultl-olo Inlet Vortex Valve 

A multiple inlet vortex valve is a cylindrical cavity with two 
opposing sets of tangential inlets arranged so that with either set 
flowing alone maximum vortex strength is obtained. Figure 1.0 
illustrates this geometry. The clockwise inlets are supplied from 
a plenum above while the counter-clockwise inlets are supplied 
from a plenum bole::. The interesting feature provided by the 
geometry is its input impedance. Input impedance is defined as 
the slope of the curve of pressure versus mass flow at the counter- 
clockwise ports while the clockwise ports are supplied with constant 
mass flow through a choked orifice. The resulting characteristic is 
as shown in Figure 1.1 Positive flow is defined as flow into the 
vortex cavity. 





FIGURE 1.0 SCHEMATIC OF A MULTIPLE INLET VORTEX VALVE 



( FLOW ) 



- 4 - 




FIGURE l.l TYPICAL IMPEDANCE CHARACTERISTIC 



A discussion of this characteristic vd.ll lead to a better 



unders tending of the operation of the valve . Suppose that the bias 
flow set of inlets cause a clockwise vortex while flowing alone and 
therefore the other inlet set would provide a counter-clockwise 
vortex wliile flowing alone. For discussion purposes, call the 
counter-cloclr.-d.se inlets number 1 and the clockwise inlets number 
2. Referring to Figure 1.1, region 1 corresponds to flow out of 
the vortex cavity through the number 1 inlets. Point "B" corre- 
sponds to the pressure caused by the flow through the number 2 
inlets with no flow thorough the number 1 inlets. As flow is 
introduced into the number 1 inlets , it decreases the vortex 
strength and causes the pressure to decrease; thus, region 2 of 
the characteristic is generated . At point "A" the angular momenta 
of the number 1 and number 2 inlets cancel. As more flow is 
introduced into the number 1 inlets , a counter-cloclad.se vortex 
is formed and the pressure increases; thus, region 3 of the 
characteristic is generated. 

This characteristic provides a negative resistance region and, 
therefore, a relaxation oscillator could be built. A flow source, 
a horizontal line on the coordinates of Figure 1.1, would have 
stable intersections in all regions of the characteristic . How- 
ever, a pressure source , a vertical line on the coordinates of 
Figure 1.1, would have stable intersections in regions 1 and 3 and 
unstable in region 2. Suppose there is a disturbance in the flow 



- 6 - 



source which causes a flow increase momontarj ly while the pressure 
rejoins fixed, in region 2 the pressure required by the number 1 
inlets decreases and, therefore, the source provides none flow 
until a stable intersection is reached vertically above in region 

3. 

A relaxation oscillator can be constructed if the nurber 2 
inlets are supplied through a chokod valve with a constant bias 
flow while the number 1 inlets are supplied through a volume. 

The flow to the volume is set by a choked valve at appr o:cu '.ate ly 
one-half the bias flow. The volume acts as a pressure source. 

This combination provides excessive flow for region 1 and 
insufficient flow for region 3* The result is a limit cycle in 
which the volume discharges in region 3 until the vortex valve 
switches to region 1 where the volume charges until the vortex 
valve switches back to region 3, and the cycle repeats. Figure 
1.2 is a schematic diagram of the relaxation oscillator. 

1.4 S cone 

The scope of this thesis is to design and construct a proto- 
type to examine one parameter of the internal geometry of the 
valve, to predict and measure the pressure sensitivity of the 
oscillator, and to predict and measure the temperature sensi* x— 
ity of the oscillator. 
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FIGURE 1.2 



SCHEMATIC OF RELAXATION OSCILLATOR 



APPARATUS 



2.1 Valve Design 



Host of the design lias been done empirically due to the 
absence of analytical methods . The valve has been designed so 
that the ratio of inlet area to outlet area could- be varied . 
Previous designs, see Reference (1), were constructed so that the 
opposing tangent! 2.1 inlet sets were lying in the same plane . 
Therefore, the present design separated the sets of tangential 
inlets axially by 0.25 inches. Both of the geometric variations 
mil provide useful empirical information for further design and 
possibly lead to a correlation for analytical design. 

The vortex chamber diameter was chosen to be 2.0 inches. 

The vortex cavity has two outlets; each outlet area, is 0.109 so. 
inches. To provide versatility in the inlet conf iguration , 
twenty-four inlet slots were provided every fifteen degrees 
around the vortex chamber perimeter . In order to provide twenty- 
four inlets aroung the circumference of a two inch diameter circle 
and have part of the circumference remaining, the inlets must r.?J:e 
an angle of greater than 35 degrees with a tangent to the circum- 
ference of the circle. The 0.125 X 0 . 3.25 inch inlets were machined 
at a 40 degree angle to a tangent at the vortex chamber perimeter . 
Figure 2.0 is a detailed drawing shewing only 4 inlets for clarity. 
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SCALE : l/20"= 1/8" 



FIGURE 2.0 DETAILED DRAWING OF A VORTEX VALVE 

PNEUMATIC OSCILLATOR 
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The inlets a re supplied iron plenums above and below part 2. 

The volume was designed as an integral part of the valve . 

This design assured equality of the volume pressure and the pressure 
at inlets number 1. If the volume is chosen very Large compared to 
the vortex cavity volume , the amplitude of oscillation corresponds 
closely to the knees of the curve shown in Figure 1.1. The reason 
for this is that the charging time of the volume is much greater 
than the switching time of the valve. As the volume approaches 
the volume of the vortex cavity, the oscillation stops; presumably 
due to the increasing influence of the valve dynamics. The present 
device was designed with the volume nearly 4 times larger than the 
volume of vortex chamber cavity to assure oscillation. The volume = 
5.85 cu. in. 

Figures 2.1 and 2.2 show a photograph of both sides of part 1 



in Figure 2.0 
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FIGURE 2.1 BIAS SIDE OF VORTEX VALVE 




FIGURE 2.2 VOLUME SIDE OF VORTEX VALVE 
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2.2. Test Apparatus 

Figure 2*3 is a schematic of the experimental setup. Figure 
Z.k- is a photograph of the apparatus . Helicoidal tost gages were 
used to measure pressure and Fischer-Porter variable area motors 
were used to measure flow rates. The air heater was a b.5!cw 
Wiegand gas heater. Temperatures were measured vrith copper- 
constantan thermocouples. Pressure oscillation was measured using 
a Eynisco PT 25"25 pressure transducer. Frequency was measured by 
stopping a Lissajous figure which was generated with the transducer 
output on the vertical amis and a sine function generator output on 
the horizontal axis of an oscilloscope. The frequency was counted 



using a Hewlett-Packard frequency counter. 
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FIGURE 2.3 EXPERIMENTAL SETUP 



0 





FIGURE 2.4 APPARATUS 
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CHAPTER III 



EXPERDENIA L DETEJC-ZMATIOil OF STATIC VALVE CHARACTERISTICS 



3 • 1 Ifcroerimental Pr o cedur e 

Referring to Figure 2.3, the static characteristics of the 
vortex valve were obtained as follows : 

a) The flow rostering pressure was set and held 
constant at UO psig. 

b) Bias flow, "C" , was set and held constant at 
either 11.1 SCFM or 5*35 3CFK. 

c) Region 1 of the characteristic curve was obtained 
by opening valve 11 3” and recording flow at 
flowmeter "S’ 1 while valve ,r A M remained closed. 

d) Regions 2 and 3 of the characteristic curve 
were obtained by opening valve "A" and record- 
ing flow at flowmeter "A" while valve B 
remained closed. 

e) The pressure at the number 1 inlets was 
measured with a mercury manometer. 

The above procedure was repeated for each of the tost conditions 
shown in Table 3*1* 

Leakage flow vxas 1.6% of a total input of 28.25 SCFM. 



3.2 Data 



‘- 16 - 



Figures 3.1 to 3*S show the static valve characteristics for 
the test conditions of Table 3*1« 



TABUS 3*1 
TEST CONDITIONS 



ist No. 


h A 


Bias Flow (SCFK) 


Air Temperature (°?) 


T 1 


0.30 


5.25 


65 


T 2 


0.30 


5.25 


250 


T 3 


0,45 


5.25 


65 


T 4 


0.45 


5.25 


250 


T 5 


0.60 


5.25 


65 


T 6 


0.60 


>.25 


250 


T 7 


0.91 


5.25 


65 


T 8 


0.91 


5.25 


250 


T 9 


0.30 


11.10 


65 


T10 


0.30 


11.10 


250 


Til 


0.45 


11.10 


65 


T12 


0.45 


11.10 


250 


T13 


0.60 


11.10 ' 


65 


T14 


0.60 


11.10 


250 


T15 


0.91 


11.10 


65 


T16 


0.91 


11.10 


250 
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FIGURE 3.S STATIC VALVE CHARACTERISTICS 
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FIGURE 3.G STATIC VALVE CHARACTERISTICS 
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FIG UR E 3.7 STATIC VALVE CHARACTERISTICS 
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FIGURE 3.8 STATIC VALVE CHARACTERISTICS 
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3 • 3 Analysts for Stati c Char acte rist ics 



1 he analysis for this do vice was originated by Ir. Ado r l C. 
Eell t Reference (l). To provide continuity to this report, it is 
repeated here . 

See Figure 3«9* The governing equations for the device are 



as follows : 



a ) Source Equation 

= constant 



b) Volume Equation (continuity) 



W. - W, 
in 



/oV dt> 
B dt 



c ) Load 



W 1 ~ g (P]_* T, goo. s try) 



( 1 ) 

( 2 ) 

(3) 



The volume is an integral part of the valve ; therefore , 

= p. Since W- ^ cannot be determined as an analytic function of 
pressure, p, the equation (2) cannot be integrated directly. 
Therefore , a graphical technique will be used for the solution . 
The technique is slope-line integration - see Appendix B. 

Rewriting the volume equation in finite difference form I 



Sp = (u in - wp 5T 



(M 



ST = J3_ St 



where 



(5) 
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FIGURE 3.9 SCHEMATIC OF RELAXATION OSCILLATOR 



This fora of the equation states' "the change in pressure over 
an increment of time is equal to the average not flow into the 
chamber during that time multiplied by an integrating tine factor" 
The time constant of the volume is given by t 

T = 

B 

Frequency of oscillation is inversely proportional to the 
time constant of the volume . Since 

B = kp 

and p _ i 

p rt 

the tiro constant of the volume is inversely proportional to 
temperature. If the static characterist5.es were independent of 
temperature , we could predict that frequency would be directly 
proportional to temperature . Since the static characteristics 
do vary with temperature , the temperature sensitivity prediction 
becomes more involved. Furthermore, the frequency is insensitive 
to pressure . However , this is only true if mass flows arc hold 
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3 *4 I frnvvr.ilc B ehav ior , The o retical 

The predicted temperature sensitivity was obtained by assum- 
ing a linear variation between 65° F and 250 c F. The frequency at 
the end points was predicted using the graphical technique 
mentioned in Section 3*3 and illustrated in Figure 3* 10. Frequency 
is calculated by relating the ST period of oscillation to real 
time through equation ( 5) • Pressure history can also be predicted 
by plotting pressure versus 

For purposes of illustrating the method, the predicted fre- 
quency of oscillation for Figure 3*10 is calculated. Figure 3*10 
shows .the period of oscillation to be 25 ST . The slope of the 
lines is: 



ST = 2.0 psig 



SCFM 



( 9) 



Relating this period of oscillation to real time 



St = V ST 
kRT 



( 10 ) 



where V = 3*3 85 x 10*" 3 c\i. ft. 

k = 1.4 

T = 525 °H 

R = 1.72 x 10 3 ft. 2/ sec. 2 op v 
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6t = 49.5 k sec. therefore 

f = 20.2 Hz 



The results of the 



integration of the remaining static character- 



istic curves are shown in Figure 3»H 
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TEST 9 
Wj n = 7. 1 7 SCFM 




FIGURE 3.10 SLOPE LINE INTEGRATION TECHNIQUE 



FREQUENCY (Hz) 



XD 




FIGURE 3.11 PREDICTED TEMPERATURE SENSITIVITY 
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CII/iPTBH IV 
EXPERH-iElffAL RESULTS 

4.1 Geometric Eff ects • 

The geometry of the valve has many effects on its static 
characteristics. The static characteristic changes considerably 
as the inlet to outlet area ratio increases. The characteristic 
also changes with temperature. As the characteristic changes, 
so does the frequency and the output waveform. See Figures 4.1 
and 4.2. Furthermore , the temperature sensitivity of the device 
can be changed by changing valve geometry as shorn in Figure 
3.10. The geometric configuration where A^/A q - 0.60 predicts 
a 34.6 percent change in frequency for a 35*2 percent change in 
absolute temperature . 

The position of the inlet slots around the vortex chamber 
were varied holding Ag/A 0 constant. It was found that spacing 
the inlets symmetrically or assvmmstrically had little effect 
on the valve characteristic , provided the two sets of tangential 
inlets remained axially opposed. 

The frequency of oscillation should be inversely proportion- 
al to the volume for a given vortex cavity geometry. Scaling down 
of the entire device should increase f requeue;/ if the laws of 
dynamic similitude apply. However, it is not clear how scaling 
down the vortex cavity effects the static characteristic of the 



device 
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Aj/A 0 = 030 
W jf / BIAS = 0.646 
6.9 - 13 4 PSIG 
TIME - 20 msec/cm 




Aj/A 0 = 045 
W in / BIAS = 0.505 
43 - 9.8 PSIG 
TIME = I Omicc/cm 



FIGURE 4.1 VOLUME PRESSURE HISTORY 
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Aj /A 0 = 0 60 
Wi,/ BIAS = 0.675 
3 3- 8 1 PSIG 
T IME = 10 msec/cm 




FIGURE 4.2 VOLUME PRESSURE HISTORY 
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4.2 Pres sure Sensit ivity 

The device is very nearly insensitive to supply pressure if 
air temperature remains constant. The actual pressure sensitivity 
of the device is shown in Figure 4.3. The sensitivity was measured 
by setting flow conditions at a metering pressure of 40 psig. 

Supply pressure was then increased to 80 psig, then decreased to 
10 psig* The ratio of flow into the volume to the bias flow 
remained nearly constant throughout the range tested. 

4.3 Temperature Sensitivity 



The ratio of predicted temperature sensitivity to actual 
temperature sensitivity was 3*2 to 1 for a given geometry. The 
main problem encountered while attempting to measure temperature 
sensitivity was' a controllable constant temperature source . The 
fact that such apparatus wa.s not available made it impossible to 
obtain steady state temperature data. A further shortcoming of 
the apparatus was that the device was constructed of brass and, 
therefore, provided a large heat sink. 

However, temperature sensitivity data was obtained ’bile the 
air temperature was steadily increasing. Though this is highly 
undesirable, it was the best information capable of being produced 
from this apparatus. To guard. against misrepresent:: ng the temp- 
erature sensitivity ; temperature was measured at the valve inlet. 
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FIGURE 4.3 PRESSURE SENSITIVITY 
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valve outlet, and in the large volume . Table 4.1 clearly shows 
the large temperature gradients present . 

Figure 4.4 shows the actual temperature sensitivity of the 
device, where temperature of the air in the volum is plotted. 

i 

Although the sensitivity does not agree with that preducted in 
Section 3*4, the author feels that the undesirable thermal 
characteristics of the valve are a major cause of deviation. 

The predicted gain was obtained with the entire device at either 
65° F or 250° F. However , steady state temperatures were not 
obtained during the temperature sensitivity tests. Further 
pursuit of temperature sensitivity data with this apparatus was 
not attempted. 
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AIH 



TABLE 

TSITKRATUi<E3 



Inlet Temo ( °lQ Ou tlet Tern (°F) 



Volr-ire T out p ( °F ) 



65 


65 




80 


77 


70 


97 


92 


80 


116 


106 


90 


126 


11? 


100 


139 


127 


no 


152 


140 


120 


162 


150 


130 


171 


159 


140 


101 


I/O 


150 


191 


181 


160 


200 


190 


170 


210 


201 


180 


219 


211 


190 


226 


220 


200 


235 


229 


210 


2 45 


240 


220 


253 


251 


230 


262 


262 


240 


270 


270 


250 



FREQUENCY (Hz) 



_o< 

-3c 



A/A^O.3 
BIAS = 11. 1 SCFM 
W in = 7.17 SCFM 
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FIGURE 4.4 TEMPERATURE SENSITIVITY 
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cHAprait v 

CONCLUSIONS AiiD SUGGESTIONS 

Geometric parameters have a significant influence on the valve 
characteristics. Herein lies an area where much more investigation 
should be conducted. The effect of each of the geometric parameters 
should be studied to lead eventually to an optimum design. The 
effect of scaling the device is vitally important to its eventual 
application and , therefore , should be studied . 

This work has shown the device to be in fact temperature 
sensitive and pressure insensitive. Future work should begin by 
constructing a mold from the existing device and casting a new 
valve using a non-heat-conducting material, possibly plaster. 

Prior to attempting to measure the temperature sensitivity , a 
controllable temperature source must be obtained. If the results 
are acceptable , a device should be constructed of a material which 
follow the guidelines established in Appendix A. If a desirable 
temperature sensitive device is designed, it must be tested for 
speed of response to step changes in temperature. 
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APPEHDTX A 

HEAT TRANSFER EFFECTS 



To bo a useful temperature sensor, it is desirable to sense 
the temperature of the air entering the device. The following is 
an order of magnitude study to determine a material which results 
in the least heat transfer to the device. 

For this study, assume a solid homogenous body whose physical 
properties ere constant; it is initially at a uniform temperature 
Tj_, and is suddenly plunged into a fluid whose temperature is 
uniform at T^. ; the heat transfer coefficient at the surface of 
the body remains constant. 

Assume that the Biot number (^/hr 0 ) is greater than 6. With 
this assumption, the internal temperature gradients are negligible , 
i. e. the plate is thin. 

An energy balance may be written for the body: 

Heat input = increase of internal energy 

d.Q = Ah (T f -T) dt = c/oV dT (A.l) 

With the assumptions made, integration of (A.l) yields: 

t = cpY In T f - T-l (A. 2) 

' hA T f - T 2 



Equation (A. 2) gives the 
temperature T^ to T 2 . 



time required to heat the body from 



Solving equation (A.2) for t’nc tine to rise 63.2 percent of 
the way between an initial and final tei,:perature yields the tine 
constant : 



t 




(A 



For the temperature sensor, we want the material to respond 
to a change in temperature quickly* 

For flowing gases, the order of magnitude for h is 2 ~50. 

The actual value vd.ll depend on the gas. The time constant decrees 
es when the device is scaled down since V/A decreases with scaling. 
The best material to use to construct the device vail have the 
smallest pc product vrith the constraint that the Biot number is 
greater than 6. Table A.l shc*^ the therrr.ophysical properties «t 
1112° F. The half thickness, r D , of the material for this type 
device is on the order of 0,02 ft. Assuming h to be 50, the 
thermal conductivity, k, for a Biot number equal to 6 is on the 
order of 6. The pc product for the various alloys listed in Table 
A.l are the same order of magnitude. 

The design of an operational device could isolate the device 
from its environment by constructing a plenum chamber around the 
entire device and discharge the exhaust gas from the valve into 
the plenum. This design should greatly reduce temperature 
gradients in the device. 
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TABIE A.l 



THEPJ'Or HYSICAL PROPERTIES AT 1112° F 



pc (BTU/ft.3 op) k (STU/lir.ft.°F) 



Copper, piare 


51.1 


20^ 


Iron, pure 


53.3 


23 


I'blybdenum 


33.2 


6l 


Nickel 


59.0 


32 


Steel, mild 


55.0 


19 


Stainless steel 


53.5 


13 


TTingsten 


38.6 


. 65 



The author feels a device constructed of stainless steel 
(k - 13) and designed to be isolated fron its environment should 
minimize heat transfer to the device and thereby permit the device 
to be used as a sensor of input air temperature* 
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APPENDIX B 

SLOHS-LINS INTEGRATION 

This appendix describes slope -lino integration in genera?, 
and explains its use in analyzing the static valve characteristics. 
Reference (6) is the source of the general description of the 
technique . 

See Figure 3, 1. Suppose we are given two rectangles, 
and A£, having equal bases Ax. The areas of the rectangles are 
therefore proportional to their heights, h_ and h,, 

fi = b (B.n 

A 2 h 2 

Suppose we are also given two inverted isosceles triangles 
each having the sane vertex angle. The bases, a-^ and of the 
inverted isosceles triangles are also proportional to their 
heights, and h£» 



^1 ~ \ (3.2) 

a 2 h 2 

Therefore, the bases of the inverted isosceles triangles are 
proportional to the areas 
A^ and A^. 



of the corresponding rectangles 






r ! 




FIGURE B.l 





FIGURE B.2 
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(B.3) 



Suppose J~ 2 _ and were superimposed as shown in Figure B.2. 
They are in effect inscribed in and circumscribed about a traps - 



The area can be found graphically; therefore , by adding half 
of the "V 11 used to determine A^ to half of the "V 11 used to deter- 
mine o The result is a "V" with a linear distance a^ which is 



A curve can be approximated by a series of trapezoids; 
therefore, the area under the curve can be found graphically. If 



all having equal vertex angles will determine the area of tho 
trapezoids# See Figure D#3* Note that the points in the slope 



y, not values of x. 

Suppose we assume that the procedure is exact. The assump- 
tion is valid if the curve is a polygon rather than a smooth 
curve. The procedure is more precise for a smooth curve as the 
width of the trapezoids approaches zero, i. c. steep slopes. 




a measure of area Ky See Figure B.2. 



the widths of the trapezoids are all the same, a series of Vs 



line diagram record the growth of area plotted against values of 
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FI6URE B.3 




FIGURE B.4 
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Since 


dA = ydx 


(3.4) 




y = dA 


(3.5) 




dx 





Theroforo , Figure B.3 can be viewed as ’ 



A' = dA = f (A) (B.6) 

dec 

Such a graph is a phase-plane plot. 

Now , let us suppose we were given a graph of A* as a function 
of A and we wish to obtain A as a function of x. V’e first choose 
a Ax interval. Vfc then can set up the two slope-lines that 
constitute the V. In Figure B«4 we see that AA = A* Ax since the 
area of a rectangle A' Ax can be represented by AA. We now start 
at any point on the curve and. construct the slope •'Hanes fror; the 
given curve to the axis and back to the curve , etc . We obtain a 
series of points on the curve spaced Ax units apart. V’e can now 
calibrate these points with appropriate values of x. When we have 
finished, we can plot A as a function of x. We have successfully 
solved the differential ecfuation (B.6). We can start anywhere 
along the curve since this reflects the arbitrary choice of the 
constant of integration. 

The direct application to the analysis for the static valve 
characteristic is discussed in tine following paragraphs. The 
analogy that can be drawn between the present analysis end the 



previous discussion is straightforward. The static valve charac- 
teristics together with the flow into the volume, Uj n , represent 
a phase-plane plot. Since the not flow rate into the volume, 

V/ in - V7-J =y dp is plotted against p. Therefore, the analogous 
dt 

quantities are • 



A 


P 


X 


t 


dA 


dp 


ebe 


dt 



Since we are given a graph of dp as a function of p, we can plot 

dt 

p as a function of t and , therefore , can obtain the frequency of 
oscillation . 

The procedure is to establish slope-lines of slope ST/2 as 
shown in Figure 3*1° and construct the zig zag pattern between 
W^_ n and The linear distance between the Vs represent Sp. 

We obtain a series of points on the curve spaced ST units apart . 
We can now plot p as a function of T and thus determine the 
period of oscillation. The period can be related to real time 
through the equation 

St = _V_ ST 
kRT 



and thus the frequency of oscillation can be computed. 
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